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ABSTRACT 


A  totally  merged  MOSFET  and  BJT  has  been  proposed.  The  device  exhibits 
complicated  non-linear  characteristics  under  certain  operating  conditions.  Due 
to  the  distributed  character  of  this  novel  merged  device,  a  straightforward 
lumped  device  approach  Is  not  adequate.  A  distributed  model  Is  proposed  and 
analyzed  using  numerical  techniques.  The  active  device  area  Is  shown  to  be 
affected  by  bias  and  contributes  to  the  non-linear  character  of  the  character¬ 
istics  under  suitable  conditions.  Several  gate  shapes  are  treated. 


1.  Introduction 

A  novel  merged  MOSFET/BJT  device  (the  MOSBJT)  was  previously  described  [1]. 

The  device  is  fully  merged  which  results  in  a  distributed  behavior.  The  MOSFET 
channel  can  be  used  as  an  emitter  or  as  a  collector.  In  order  to  analyze  the 
device,  a  distributed  FET  BJT  model  must  be  used.  Approximate  analytical  models 
have  been  treated  [2]  and  will  be  reported  elsewhere  and  is  similar  to  that 
described  in  Report  III. 

The  MOSBJT  numerical  model  is  based  upon  an  equivalent  circuit  for  the  MOSBJT. 
The  construction  of  the  equivalent  circuit  is  as  follows:  First  the  distributed 
bipolar  injection  into  the  MOSBJT  channel  is  modeled  by  dividing  the  device 
length-wise  into  many  equal  sections  (refer  to  Fig.  3.1a,  b).  Each  section  is 
then  represented  by  a  n-channel  metal  oxide  semiconductor  field  effect  transistor 
(MOSFET)  connected  to  a  npn  bipolar  junction  transistor  (BJT)  resulting  In  the 
equivalent  circuit  shown  in  Fig.  3.1c.  Built  into  this  circuit  model  is  the 
assumption  that  minority  carrier  transport  occurs  vertically. 

The  npn  BJTs  and  the  n-channel  MOSFETs  in  the  MOSBJT  equivalent  circuit 
(Fig.  3.1c)  are  represented  by  their  classical  DC  models.  These  are  the  Ebers-Moll 
BJT  model  [3]  and  the  distributed  analysis  MOSFET  model  [4], 

The  Ebers-Moll  BJT  model  represents  the  npn  BJT  with  the  equivalent  circuit 
shown  in  Fig.  3.2.  From  the  equivalent  circuit  the  following  expressions  for  the 
collector  current  (I£)  and  the  base  current  (Ig)  for  the  transistor  are  determined. 


qVRF/KT  qVRf,/KT  qVRr/KT 
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(3.1) 

(3.2) 


The  Ebers-Moll  model  parameters  (ap.aR,  leg,  Ics )  for  the  BJT  elements  in 
the  MOSBJT  circuit  model  are  extracted  from  tne  experimental  data  of  a  typical 
rectangular  MOSBJT  operating  with  the  entire  channel  active.  The  magnitude  of 
the  BJT  collector  current  for  a  particular  section  (lex)  proportional  to  the 
collector  area  of  that  section.  The  proportionality  constant  is  determined  by 
normalizing  the  area  of  each  section  to  the  total  collector  (channel)  area.  This 
factor  Is  then  used  to  scale  the  total  collector  current  expressed  by  Eq.  (3.1) 
yielding: 


I 


Cx 


Vc 


(3.3) 


where 

K  _  collector  area  of  the  xth  section 
\  ~  total  collector  area 

x  =  Integer  specifying  the  section  of  the  numerical  MOSBJT  model  being 
referred  to 

Similarly  an  expression  for  the  base  current  contribution  of  a  particular  section, 
x  can  be  expressed  by  Eq.  (3.4). 


2 


(a)  Top  view 


Emitter 


(b)  Cross-sectional  view 


Emitter- 


Fig.  3.1  Development  of  the  fiOSBJT  equivalent  circuit,  (a)  Top 
view  of  a  rectangular  MOSBJT.  (b)  Cross-sectional  view 
of  a  rectangular  MOSBJT.  (c)  Equivalent  circuit  of  the 


(3.4) 


where  Kx  and  x  are  defined  as  above. 

The  n-channel  MOSFETs  are  modeled  with  the  distributed  MOSFET  model  and  the 
gradual  channel  approximation  [3].  Incorporated  into  the  MOSFET  model  are  the 
assumptions  that  the  channel  voltage  can  be  approximated  by  a  half  of  the  sum  of 
the  source  and  drain  voltages  (^(Vs  +  Vq))  and  that  the  channel  length  is  much 
greater  than  the  depletion  region  at  the  drain.  The  MOSFET  model  yields  the  fol¬ 
lowing  expression  for  the  drain  current  (Id). 

!0  =  yn  L  Cox^VG  “  VT  "  (  °2  S^^VD“VS^  ^3'5^ 


where 

W  =  width  of  the  channel 
L  5  length  of  the  channel 
yn  =  electron  mobility  in  the  channel 
Vj  i  threshold  voltage  . 

Vq  =  applied  drain  voltage 
V$  i  applied  source  voltage 
Vq  =  applied  gate  voltage 

Solving  Eq.  (3.5)  for  the  drain  voltage  (Vq)  yields  Eq.  (3.6). 
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The  n-channel  MOSFET  elements  in  the  MOSBJT  circuit  model  (Fig. 
represented  by  Eqs.  (3.5)  and  (3.6).  Since  each  MOSFET  represents  a 
the  MOSBJT  channel,  the  length  of  each  MOSFET  should  be  equal  to  the 
channel  length  divided  by  the  number  of  sections  as  expressed  by  Eq. 


3.1c)  are 
section  of 
total  MOSBJT 


(3.7). 


L’  =  L/N 


(3.7) 


where 

L‘  =  the  channel  length  for  the  MOSFET  circuit  elements 
N  =  the  number  of  sections  in  the  MOSBJT  circuit  model 

For  a  rectangular  gate  MOSBJT  the  width  of  the  gate  of  the  MOSFET  circuit 

element  will  simply  be  equal  to  the  width  of  the  MOSBJT  gate.  MOSBJT’s  with -  - 

non-rectangular  gate  geometries  are  easily  simulated  by  representing  each  section 
with  a  rectangular  MOSFET  of  length  L/N  and  width  equal  to  the  average  width  of 
the  MOSBJT  channel  in  that  section. 

The  computer  algorithm,  listed  in  Appendix  A  and  outlined  in  Fig.  3.3,  consists 
of  4  modules.  The  function  of  these  modules,  labeled  A  through  D,  are  as  follows: 

In  Module  A,  the  MOSBJT  material  properties,  device  characteristics,  geometry,  and 
operating  bias  are  defined  and  various  required  parameters  calculated.  The  purpose 
of  module  B  is  to  determine  the  extent  of  channel  cut-off.  Next  in  module  C  the 
spatial  dependence  of  voltage  and  current  along  the  channel  is  determined  and 
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Fig.  3.3  Flow  chart  of  the  numerical  MOSBJT  model 


recorded.  Finally  in  module  0  the  base  current  (Ig)  is  calculated  and  its  value 
along  with  the  rest  of  the  data  is  printed  or  plotted. 


In  module  A,  the  numerical  MOSBJT  model  parameters  are  defined.  These  parameters 
are  listed  in  Table  3.1  along  with  their  typical  values.  Values  for  <j>p,  C,L,  Vpg,  Vt 
are  determined  by  evaluating  the  following  expressions  derived  from  conventional  MOS 
theory  [3]. 


KT 


V  q  ln(W 
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VT  1  Vfb  +  2l+pl  +Cj7^sqNa(2|*p|-VBE)  (3.11) 

A  substantial  portion  of  the  iterative  algorithm  of  the  numerical  MOSBJT 
model  is  contained  in  modules  B  and  C.  In  these  modules  the  spatial  dependence 
of  the  MOSBJT' s  channel  voltage  and  current  are  determined  for  a  given  device 
bias  (Vpp,  Vgc,  and  Vq).  Once  the  channel  voltage  and  current  profiles  are  known 
the  terminal  drain  current  (Ip)  and  base  current  (Ib)  can  be  determined. 

The  basic  MOSBJT  circuit  model  is  shown  again  in  Fig.  3.4  with  the  various 
nodes,  node  voltages,  and  currents  labeled  with  the  following  notation: 

X 

I  MX 
N 

A  simple  approach  to  determine  the  node  voltages  and  currents  would  be  to 
start  at  the  drain  where  the  voltage  is  known  and  then  calculate  the  node  voltage 

and  current  moving  from  node  to  node  in  the  direction  away  from  the  drain.  This 

approach  does  not  yield  a  solution  because  finding  Vr  (the  voltage  at  the  node 
adjacent  to  the  drain)  requires  that  1^  be  known.  Im^  is  equal  to  the  sum  of  the 
BJT  collector  currents  at  all  other  nodes.  These  collector  currents  are  in  turn 
dependent  on  their  respective  node  voltages.  The  interdependence  of  node  voltages 
and  currents  lead  to  the  development  of  an  iterative  approach  for  the  numerical 
MOSBJT  model . 

In  this  approach  an  initial  guess  is  first  made  of  the  voltage  (V-j)  at  the 
node  furthest  from  the  drain  (node  1).  Next,  Iq,  can  be  solved  for  using  the 

expression  for  the  collector  current  of  a  BJT  obtained  from  the  Ebers-Moll  BJT 

model  (Eq.  (3.3)).  The  current  flowing  through  the  MOSFET  interconnecting  nodes 
1  and  2,  (Im)  is  equal  to  Iqi .  A  general  expression  relating  the  node  currents, 

(Eq.  (3.12))  is  derived  with  Kirchoff's  current  law. 


=  integer  indicating  the  X**1  section  of  the  numerical  MOSBJT  model 
=  current  flowing  through  the  MOSFET  located  in  the  Xth  section 
=  collector  current  for  the  BJT  located  in  the  Xth  section 
=  voltage  at  the  Xth  node 

=  the  number  of  sections  for  the  numerical  MOSBJT  model 
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Table  3.1  Typical  MOSBJT  Device  Parameters 


Parameter 

Nomenclature 

Typical  Value 

MOSBJT  channel  length 

L 

3.68E-3  M 

MOSBJT  channel  width 

W 

2.96E-3  M 

Plumber  of  sections  in  model 

N 

100 

Gate  voltage  w.r.t.  the  emitter 

VG 

30V 

Drain  voltage  w.r.t.  the  emitter 

VDE 

20  V 

Base  voltage  w.r.t.  the  emitter 

vBe 

0.38  V 

Metal  work  function  (Al) 

^m 

4. lev 

Semiconductor  work  function  (Si) 

(j>S 

4.9ev 

Oxide  permittivity 

£OX 

3.45E-1 1  F/M 

Semiconductor  permittivity  (Si) 

es 

1.04E-10  F/M 

Semiconductor  doping 

Na 

1.5E16  cm-3 

Temperature 

T 

300°K 

Oxide  thickness 

Tox 

1 .2E-7  M 

Electric  potential  of  the  semiconductor 

4p 

-0.358V. 

Oxide  capacitance  per  unit  area 

Cox 

2.88E-4  F 

Flatband  voltage 

VFB 

-1.08V 

Threshold  voltage 

vT 

1.72V 

Ebers-Moll  model  parameter 

aF 

0.9615 

Ebers-Moll  model  parameter 

aR 

0.9975 

Ebers-Moll  model  parameter 

rES 

1.04E-8 

Ebers-Moll  model  parameter 

!cs 

1.002E-8 

-i 


ft 


%  \  %  % 
k-i-A  5.  i*i A  ■  «4 


!mx  =  !MX-1  +  1 


(3.12) 


CX 

Knowing  V-j  and  ty,  V2  is  then  calculated,  using  Eq.  (3.6)  in  which  a  MOSFET's 
drain  voltage  (Vx+i)  is  expressed  in  terms  of  the  MOSFET's  source  voltage  (Vx)  and 
drain  current  (ImxI*  This  process  continues  until  the  voltage  at  the  drain  (Vdcal) 
is  found.  If  Vqcal  1S  not  equal  to  the  applied  drain  voltage  (Vde)  iterative 
adjustments  are  made  to  the  initial  guess  until  the  desired  accuracy  is  achieved. 

From  the  basic  understanding  of  the  MOSBJT's  operation  [Final  Report  I] 
it  is  known  that  under  certain  operating  conditions  a  portion  of  the  MOSBJT  channel 
may  be  "cut-off"  -  incapable  of  collecting  the  injected  base  minority  carriers.  To 
incorporate  channel  "cut-off"  into  the  numerical  MOSBJT  model,  the  iterative 
algorithm  is  partitioned  into  two  modules  B  and  C. 

In  module  B  the  extent  of  channel  "cut-off"  is  determined.  The  algorithm 
contained  in  module  B  is  diagramed  in  Fig.  3.5.  First  the  minimum  channel  voltage 
with  respect  to  the  emitter  (VcHEMIN) >  at  which  the  channel  is  still  active  as  a 
collector  is  found.  This  is  done  by  first  solving  for  the  collector-base  voltage 
( vBc)  at  whjch  channel  "cut-off"  occurs.  Since  channel  "cut-off"  occurs  when  the 
BJT  collector  current  (lex)  1S  equal  to  zero,  Eq.  (3.3),  the  expression  for  lex 
can  be  set  equal  to  zero,  as  shown  in  Eq.  (3.13),  and  then  Vbc  solved  for. 

qVRF/KT  qV'/KT 

0  =  Ia  =  Kx  FIES(e  BE  -1)  -  Ics(e  BC  -1)  (3.13) 

The  sum  V^c  +  Vbf  is  the  channel  voltage  with  respect  to  the  emitter  necessary  for 
channel  "cut-off11  to  occur.  Therefore  VcHEMIN  can  be  approximated  by  adding  a 
small  voltage  to  this  sum  as  shown  in  Eq.  (3.14). 


VCHEMIN  =  VBC  +  VBE  +  0,01  V 

After  Vchemin  is  found  it  is  applied  to  the  node  located  furthest  from  the  drain. 
Equation  (3.3),  (3.6),  and  (3.12)  are  used  to  calculate  the  node  voltages  and 
currents  along  the  channel  until  the  drain  Is  reached  where  VdcaI,  a  Vr+i •  Then 
Vocal  is  compared  with  the  applied  drain  voltage  Vde*  If  VocAl,  ls  greater  than 
Vqe  the  implication  is  that  the  initial  guess  is  too  large.  Since  the  voltage 
applied  to  node  1  was  VcHEMIN  (the  minimum  voltage  for  which  a  node  will  be  "active"), 
node  1  must  be  biased  into  "cut-off."  The  "cut-off"  node  is  modeled  by  setting  the 
node  voltage  (Vi)  to  (Vbc+Vbe)  and  the  BJT  collector  current  (Ici)  and  the  MOSFET 
current  (Imi)  equal  to  zero.  Then  VcHEMIN  1s  applied  to  the  next  node  (node  2)  and 
the  process  is  repeated.  This  continues  until  VOCAL  is  less  than  Vqe  upon  which 
module  C  is  executed. 

In  module  C  the  voltage  at  the  "active"  node  furthest  from  the  drain  is 
adjusted  until  VocAL  is  approximately  equal  to  Vde  to  within  a  predetermined 
accuracy.  When  this  occurs  the  MOSBJT's  terminal  drain  current  (Id)»  given  by 
IMN  and  the  MOSBJT's  channel  voltage  and  current  as  a  function  of  position  along 
the  channel  are  recorded.  Next  module  D  is  executed  and  the  MOSBJT's  base  current 
(I[j)  is  calculated.  Ib  is  found  by  summing  the  base  current  contributions  of  each 
section  as  shown  below. 
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!b  = 


N 

I 

X=1 


BX 


(3.1 


where 


N  =  the  total  number  of  sections  in  the  numerical  model 
IRv  =  the  base  current  contribution  of  the  X™  section.  An  expression  for 
IBX  is  Eq.  (3.4). 

The  numerical  MOSBJT  model  was  used  to  simulate  a  rectangular  MOSBJT  with 
the  device  characteristics  listed  In  Table  3.1.  Figures  3.6,  3.7,  and  3.8  are 
plots  of  the  theoretical  spatial  dependence  of  channel  voltage  (Vx)  and  current 
(Imx)»  and  the  base  current  contribution  of  each  section  (I|jx)»  respectively. 

In  Fig.  3.8  it  Is  evident  that  there  exists  a  slightly  greater  base  current 
contribution  from  the  sections  located  In  the  "cut-off"  region.  In  the  "cut-off" 
region  both  the  emitter-base  and  collector-base  junctions  are  forward  biased. 

This  results  in  an  increase  in  the  base  minority  carrier  concentration  and  there¬ 
fore  an  increase  in  the  amount  of  base  minority  carrier  recombination  which  is 
a  component  of  base  current.  The  forward  biased  collector-base  junction  of  the 
cut-off  region  has  an  additional  base  current  component  resulting  from  the  "reverse1 
injection  of  minority  carriers  from  the  base  into  the  channel,  recombination  at  the 
MOS  surface  and  in  the  FET  depletion  region. 

For  the  simulated  rectangular  MOSBJT  many  operating  loci  [Id(Vqe»VG»Ib)]  were 
calculated  for  a  constant  gate  voltage.  The  data  was  then  plotted  to  produce 
the  BJT-Ilke  family  of  curves  shown  in  Fig.  3.9. 

In  summary  the  numerical  MOSBJT  model  has  the  potential  of  yielding  an 
accurate  model  of  the  MOSBJT.  Although  this  model  requires  a  moderate  amount  of 
computer  computation,  non-rectangular  devices  are  easily  simulated.  The  current 
level  in  the  development  of  the  MOSBJT  is  focused  on  a  closer  examination  of  the 
MOSBJT  non-linear  electrical  properties. 

Extensive  additional  theoretical  work,  both  analytical  and  numerical  has  been 
performed.  A  portion  of  the  analytical  work  is  quite  similar  to  that  described  in 
Report  III  and  therefore  will  not  be  included  here.  It  is  anticipated  that  this 
theoretical  work  will  be  published  at  a  later  date  and  also  included  in  a  Ph.D. 
Thesis. 
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J.7  Spatial  dependence  of  channel  current  (tyx)  for  a  rectangular  MOSBJT  with  devi 
parameters  listed  in  Table  3.1. 


CURRENT  COMPONENT  OF  EACH  SECTION 


DRAIN  EDGE  PERIPHERY 

Fig.  3.8  Spatial  dependence  of  the  sectional  base  current  contribution  for  a  rectangular 
HOSBJT  with  device  parameters  listed  in  Table  3.1. 


ID-VDE  CHARACTER ISTICS  OF  THE  RECTRNGULAR  MOSBJT 


DRAIN-EMITTER  VOLTAGE  ( VDE)  volts 

Fig.  3.9  Theoretical  In-VnE  characteristics  for  a  rectangular  MOSBJT  with  device  parameters 
listed  in  Table  3.1 . 
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V  * 

*  -4 

*  THE  NUMERICAL  MOSBJT*  MODEL 

•4  -4 

*■  Developed  by  I  David  Okada  *- 

■*  Research  on  the  MOSBJT  is  being  directed  by  •* 

*  Dr.  James  Holm-Kennedv  *■ 

*  University  of  Hawaii  * 

*  Program  developed  on  a  HP  9836  with  Basic  2.0  *• 

4  -4 

*  ‘  INVENTED  BY  DR.  JAMES  HOLM-KENNDY  * 

*4  -4 


* 44# 4 444444 4-4444 444444444444444444 ** 44444444  *■* 4 ****** **** 


Main  Program 


OPTION  BASE  0 

DIM  Ic(200) ,Im(200) ,Vd(200) ,Vs( 100) , Ibi ( 200 ) 

INTEGER  N 

!  THE  MOSBJT  OPERATING  BIAS  MUST  BE  SPECIFIED 

i 

!  Vg  ...  THE  GATE  VOLTAGE 

!  Vde  . .  THE  APPLIED  DRAIN  VOLTAGE 

!  Vbe  . .  THE  EMITTER-BASE  JUNCTION  BIAS 

!  N  ...  THE  NUMBER  OF  SECTIONS  IN  THE  MODEL 

Va*20 

Vde-1  0 

Vbe -.40 

N-100 

CALL  I_v_char ( Vq, Vde , Vbe .N. Ic( *) , Im( *) ♦ Vd ( *) , Ibi ( * ) , Ibcal ) 

i 

I  THE  SUBROUTINE  I_v_char  RETURNS 

i 

!  Vd(-*)  ..  THE  ARRAY  CONTAINING  THE  NODE  VOLTAGES 

!  Ic(*)  ..  THE  ARRAY  CONTAINING  THE  BJT  COLLECTOR  CURRENTS 

!  Im ( * )  ..  THE  ARRAY  CONTAINING  THE  MOSFET  DRAIN  CURRENTS 

!  Ibi ( *) . .  THE  ARRAY  CONTAINING  THE  BASE  CURRENT  FOR  EACH  SECTION 

!  Ibcal  ..  THE  TOTAL  TERMINAL  BASE  CURRENT 

!  Im ( N )  ..  THE  TERMINAL  DRAIN  CURRENT  FOR  THE  MOSBJT 

!  Vd(N+t ) .  THE  CALCULATED  TERMINAL  DRAIN  VOLTAGE  ( VDCAL ) 

i 

FOR  L= 1  TO  N+1 

PRINT  Vd ( L ) . Ic ( L ) , Im ( L ) 

NEXT  L 
END 


-*  *  -f *  *  *  ■*  *  *  •*  *  ■*  *  *  *  *  *  *  4  «■  *  *  *  *  *  ■*  *  *  *  *#  #  #  4  ■*  *  *  *  *  *  -*  *  *  ■*  ■*  4  •#*  *  4  ■*  -*  *  *  *  -*  *  •*  *  *  *  •*  ■*  4 

Subprograms 

I_V_char  ....  Given  the  operating  bias  of  the  MOSBJT 
this  program  calculates  the  resultant 
current  and  3nd  voltage  profile  along 
the  channel  and  the  terminal  base  ana 


drain  currents. 


590 
600 

6 ! 0  !  Wid  ...  This  subprogram  to  calculate  the  average 

620  !  width  of  a  MOSFET  circuit  element. 

630  ! 

640  !  Inca  .  This  subprogram  calculates  the  ratio 

650  !  of  the  collector  area  of  a  particular 

660  !  section  to  the  total  MOSBJT  channel  area 

670  ! 

G80  1  4******#**********4**********#******-*-*-*-*-*-***-<-***-*****-*-*********-*-«-*** 

690  ! 

700  ! 

710  SUB  I_v_char(Va,Vde,Vbe. INTEGER  N, REAL  Ic(+) . Im(«) ,Vd(*> . Ibi (*) , Ibcal ) 

720  DIM  Vs(2G0) 

730  INTEGER  Nmin.Nmax.L 
740  ‘  ! 

750  GOSUB  Modules 

760  GOSUB  Moduleb 

770  GOSUB  Modulec 

780  GOSUB  Moduled 

790  GOTO  Fini 


800  ! 

gin  i  ii  i  i  i  i  i  i  i  i  i  i  i  i  i  ii  m  i  i  i  i  i  i  i  i  ii  i  i  m  i  m  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  m  i  i  i  i  i  i  i  m  i  ii  i  i  i  i  i 

820  ! 

MODULE  A 

830  ! 

840  ! 

1  . 

Define  device  parameters 

850  ! 

2. 

Define  device  geometry 

*  860  ! 

3. 

Define  material  constants 

870  ! 

4. 

Define  device  bias 

380  ! 

5. 

Calculate  necessary  intermediate  parameters 

v  190  ! 

iOO  !  !  !  !  ! 

mi . mi . . 

i  l  l  i  l  i  i  l  l  i  l  i  i  i  i  i  i  l  i  i  i  i  i  i  i  i  i  i  i  i  . . .  i  l  i  i  i  i  l 

i 

910  !  . 

920  Modulea:  ! 

930  ! 

940  ! 

Define  device 

parameters  .... 

950  ! 

960  ! 

Om  - 

METAL  WORK  FUNCTION 

970  ! 

Os  - 

SILICON  WORK  FUNCTION 

980  ! 

Eox 

-  PERMETIVITY  OF  SI02 

990  ! 

Es 

-  PERMETIVITY  OF  Si 

1000  ! 

Tox 

-  OXIDE  THICKNESS 

1010  l 

Qss 

-  SURFACE  CHARGE  DENSITY 

1020  ! 

T 

-  TEMPERATURE 

1030  ! 

Na  - 

SUBSTRATE  DOPING 

1040  ! 

1050  ! 

1060 

Flg-0 

1070 

Om-4 . 1 0 

!  V 

1080 

Os-4.9 

!  V 

1090 

Eox-3. 45E-1 1 

!  F/H 

„  1100 

Es-I .04E-10 

!  F/M 

1110 

Tox-1 . 2E-7 

!  METERS 

1120 

T-300 

!  DEG  K 

,  1130 

Na- 1 . 5E+1 6 

!  CM  *  -3 

1140 

Qss-8.E-5 

!  C/M' 2 

1150  ! 

1160  ! 

1170  !  Define  device  oeometry  ... 

11  SO  ! 


V*  "J* s'.  • 

■.■sv-.v 

o**  •.*  *;  % 

y 

* 

;i 

.  i 

•i 

j 

91 

S  1190 

1 

Lenoth-  LENGTH  OF  GATE 

v  1200 

j 

Width  -  WIDTH  OF  GATE 

■: 

^  1210 

! 

N  -  NUMBER  OF  SECTIONS 

il 

•J 

£  1220 

i 

.-i 

£  1230 

! 

,  j 

I  1240 

! 

i 

*  1 250 

Lenath  =  3 . 68E-3  !  METERS 

;?  1260 

Width=2.96E-3  !  METERS 

, 

>  1270 

i 

\  1280 

i 

■ 

til  1290 

i 

|  1 300 

j 

Specify  BJT  parameters  ... 

"  1310 

j 

‘j.  1320 

j 

AF  -  FORWARD  ALPHA 

u 

'A-  1330 

i 

AR  -  REVERSE  ALPHA 

* 

Si  1340 

j 

IES  - 

£  1350 

i 

ICS  - 

1  1360 

! 

XB  -  BASE  WIDTH 

i 

"  1370 

i 

Wj 

1380 

Af = .962 

' 

1390 

i 

Ar  = . 998 

.* 

1400 

Ies- 1 . 04E-8 

'• 

1410 

Ics-1 .002E-8 

> 

£  1420 

i 

B 

"  1430 

i 

m 

!•:  1440 

i 

I- 

•*j  1450 

i 

CALCULATE  INTERMEDIATE  CONSTANTS 

■ 

;-J  1460 

i 

1470 

i 

Vfb  -  PLATBAND  VOLTAGE 

•  ■]' 

ti  1480 

j 

Vt  -  THRESHOLD  VOLTAGE 

8 

jfl 

•  1490 

j 

D1  -  SECTIONAL  MOSFET  LENGTH 

n 

■•1  1500 

i 

Cox  -  OXIDE  CAPACITANCE/AREA 

-  /; 

.?  1510 

! 

Mu  -  MOBILITY 

1520 

i 

Op  -  ENERGY  DIFFERENCE  BETWEEN 

THE 

1530 

j 

INTRINSIC  AND  EXTRINSIC  FERMI 

ji  1540 

i 

LEVELS 

5  1550 

i 

?! 

•’  1560 

i 

;-3 

:•  1570 

I 

d 

;  1580 

Cox-Eox/Tox 

sj 

1590 

Vf b-Om-Os-Qss/Cox 

■Ij 

-  1600 

Kb- 1 .38E-23 

i  1610 

0-1 .6E-19 

1 

•:  1620 

Mu- . 05 

* 

>  1630 

N i - 1 . 45E+1 0 

1640 

Op-  - 1  * (Kb*T/G)-»LOG(  Na/Ni ) 

C 

fi]  1650 

Dl-Lenath/M 

I1:, 

U  1660 

Vt-Vfb+2*ABS(0p)+(  1  /Cox  HSQR(2*Es*Q*Na*  1 . 0E+6*(2*ABS(Op) ) )  \ 

a  1670 

i 

1 

'V  1680 

! 

IF  THE  MOSBJT  IS  BIASED  INTO  SATURATION  VDE>VG-VT  THEN  VDE 

IS  SET 

■  1690 

i 

TO  VG-VT 

V 

1700 

1 

1  , 

'•* 

'  1710 

IF  Vde>Vg-Vt  THEN 

1  . 

•  «, 

fi  1720 

Vde-Va-Vt 

|J  1730 

PRINT '"DEVICE  SATURATED  VDSAT-VG-VT  -  ",Vde 

'  | 

5  1740 

END  IF 

:•  1 750 

i 

<  1 760 

j 

-i  I77Q 

i 

PRINT  THE  ABOVE  DATA 

i  1 780 

i 

■ 

1 

m 

4 

vr*«r  ,•  *.» 

*  /  *  .  •* 

*  «,  f  'm  m  *  •  “  •  *  *.*  *.  .  ’  S  .  *  *  *  •  ’  .•  .  '  ‘  .  »  •  * 

/•  ;  •• „ 

3  1800 
•:  1310 

:•  .  1820 
1330 
1840 
x  1850 
I860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 

*1950 
1960 
v  1970 
N  1980 
1990 

;•  2000 
i  2010 


PRINT  OM= " , Om 
PRINT  "05  = " .Os 
PRINT  "EOX= " . Eox 
PRINT  "ES=" .Es 
PRINT  "T0X=" . Tox 
PRINT  "QSS=" ,Qss 
PRINT  " T= " , T , "K" 

PRINT  "NA  «" ,Na 

PRINT  "LENGTH= " , Length 

PRINT  "WIDTH-", Width 

PRINT  "NUMBER  OF  SECTIONS  IS  \N 

PRINT  " IES  =  "Jes 

PRINT  "ICS  =  ".Ics 

PRINT  "AF  =  " . A f 

PRINT  "OP  =  ",0p 

PRINT  "VFB  =  ".Vfb 


RETURN 


R  2030 

i 

2040 

! 

«  2050 

£  *2060 

j 

i 

i  2070 

i 

■  2080 

i  i  i  ii  i  i 

fr  2090 

i 

*3  2100 

Moduleb  *. 

p  2110 

j 

2120 

! 

k  2121 

j 

p  2130 

! 

P  2140 

i 

k  2150 

i 

P  2151 
K  2152 

i 

i 

h  2160 

i 

M  2170 

i 

~  2180 

i 

2190 

j 

k;  2200 

p  2210 
6  2220 
wj  2230 
“  2240 

K  2250 

i 

!!!!!!!!!!!!! m ;; i !! i; !!!!!!!!! j !!!!!! i ;; !!  ;  !!!!!!!!!!!!!!!  j  !!!!!  ; 

MODULE  B 

Determine  extent  of  channel  cut-off 

l ii it i m i i i i  i i i i i i i m i i m i i ii  ii i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  m  i  i  i  m  m  i  m  i  m  i  i  i 


Vbc  ..  The  channel  to  emitter  bias  such  that  lc*0 

Vchemim  .  The  channel  to  emitter  bias  such  that  the 
BJT  connected  to  that  node  Ic  is  slightly 
greater  than  0.  (Node  is  biased  at  the 
onset  of  cutoff) 

Leo  ...  Indicates  the  current  section  being 

examined  to  determine  if  it  is  "active" 


Vbc*-(  Kb*T/Q)-*L0G(  Af*Ies/Ics) 

Vchemin-Vbc+.OI 

Us( 1 ) -Uchemin 

Nmin*1 

Nmax-N 


Next_try:  ! 

Leo* ( Nmi n+Nmax ) /2 
Ic ( Leo ) *0 
Im  <  Leo- 1 ) *0 
Vd(Lco)*Vchemin 

i 


Loop  1 : 


FOP  L-Lco  TO  M 

IF  Q-*  ( -  Vd( L )  +Vbe )  /( Kb*T  X  -  1 00  THEN  GOTO  Around  1 


5  2360  Ic  ( L )  =  FN Inca ( Lenath ,  W  id t h  ,  N . L )  *■  (  At  *■  Ies*EXP ( Q*Vbe/  ( Kb* T  )  i  -  Ics~  ( EX?  ( 

•  (Q/(kb*T ) )M Vbe-Vd(L  )))-!)) 

^  2370  GOTO  Around2 

-  2380  ! 

■  2390  Around ! :  ! 

I  2400  Ic(L)=FNInca(Lenath,Width,N,L)-<'Af'<-Ies*EXP(Q*Vbe/(Kb*T  )  ) 

:  2410  ! 

■  2420  Around2:  ! 

2430  Im ( L ) =Im ( L- I ) +Ic ( L ) 

2440  IF  (Vg-Vt+Vd(L)  )  '2- ( 2+ ( Ug-Vt ) «Vd ( L ) +0d ( L )  *2+(2*Im<LHDl)/(Mu«FNWid 

t  (Width, Lenath. N.L  )*Cox  ) )  <0  THEN  GOTO  Cut_otf 

I  2450  Vd(L+l  )-Vg-Vt+Vd(L)-SQR(  (Vg-Vt+Vd(L) )  '2- (2*  ( Vg-Vt )  *Vd(  L )  +Vd(  L )  *2+( 

2*Im( L  ) *D i )  /  ( Mu*FNW id ( Width ,  Lenath  .  N , L  ) ■‘•Cox  )  )  ) 
j  2460  IF  Vd(L+D>Vde  THEN  GOTO  Cut_otf 

2470  NEXT  L 

<  2480  ! 

2490  ! 

|  2500  Nmax=Lco 

2510  IF  Nm i n=Nmax- 1  THEN 

•  2520  FOR  L-1  TO  Lco-1 

!  2530  Ic(L)=0 

\  2540  Im(L-H=0 

i  2550  Vd ( L ) =Gchemin 

|  2560  NEXT  L 

2570  GOTO  Found_it 

:  2580  END  IF 

|  2590  GOTO  Next_try 

|  2600  Cut_off:  ! 

2610  Nmin=Lco 

2620  IF  (Lco>-N)  THEN 

‘  2630  DISP  "  ERROR  !!!!!!  ENTIRE  CHANNEL  IS  CUT  OFF  ! H 

2640  DISP  "  RERUN  PROGRAM  WITH  A  LARGER  N  " 

!  2650  BEEP 

•  2660  PAUSE 

.  2670  END  IF 

|  2680  GOTO  Next_try 

1  2690  ! 

2700  Found,  it:  ! 

2710  '  RETURN 

2720  ! 

2730  ! 

2740  m  m  i  m  i  i  i  m  i  i  ii  m  i  i  i  ii  ii  i  i  ii  m  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  ii  i  i  i  ii  i  i  i  ii  m  ii  i  m  i  i  i  i 

'■  2750  !  . ! 

!  2760  !  MODULE  C  ! 

2770  !  j 

2780  !  !.  Iterate  until  Vdcal  =  Vde  ! 

2790  !  2.  Record  the  spatial  dependence  ot  channel  ! 

1  2800  !  voltage  and  current  ! 

2810  !  '  ! 

2820  i  i  i  i  i  i  i  i  i  i  i  i  i  i  ii  i  ii  i  i  m  i  ii  i  ii  i  ii  i  i  i  i  ii  i  i  i  i  i  i  it  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 

2830  !  ‘ 

2840  Modulec.*  ! 

2850  K-1 

2860  Loop2:  ! 

2870  FOR  L=*Lco  TO  N 

2980  IF  0-*(Vbe-Vd(L))/(Kb-»TX-l00  THEN  GOTO  AroundS 

2890  IcCL)  "FNInca  ( Lenath.  Width,  N ,  L )  *•  (  At  *Ies*EXP  (  Q*Vbe/  ( Kb*T )  )  -Ics* (  EXP  ( 

Q* ( -Vd( L ) +Vbe  ) /  ( <b*T  )  )-I  )  ) 

2900  GOTO  Around4 

29i0  Arouna3:  ! 


2920 

Ic(L) -FNInca (  Lenqth .  Width ,  N . L ) -* (  Af *Ies-*EX 

P(Q*Vbe/(Kb*T)  )  ) 

2930 

Around4:  ! 

2940 

ImlL ) =Im( L- 1 ) +Ic(L ) 

29S0 

IF  (Va-Vt+Vd(L) > '2-(2«(Va-Vt >*Vd(L)+^d<L) 

*  2+  ( 2-*Im  (  L )  *D  i )  /  ( Hu^FNW  i  d 

( H  id  th ,  Lena th  , N ,  L )  ■•’■Cox  )  )  <0  THEM  GOTO  Decrease 

2960 

Vd  C  L  + 1 )=Vg-Vt.+Vd(L)-SGR( (Vq-Vt+Vd(L ) ) '2-1 

2*(Vg-Vt)*Vd(L)+Vd(L)  ’2+( 

2*-Im(L)*Dl)/(Mu*FNWid(Hidth,Lenath,N,L)*Cox)  )  ) 

2970 

IF  Vd(L+l)>Vde  THEN  GOTO  Decrease 

2980 

NEXT  L 

2990 

Adjust:  ! 

3000 

| 

3010 

j 

2020 

IF  K=  S  THEN 

3030 

Vsm i n=Vs ( K ) 

3040 

Vsmax=Vs ( K ) 

3050 

END  IF 

3060 

IF  Fla-1  THEN  GOTO  Increase 

3070 

IF  VdT L+1 )-Vde>0  THEN 

3080 

Fla-1 

3090 

GOTO  Decrease 

3100 

END  IF 

3110 

Vsmi n=Vs ( K ) 

3120 

Vs ( K+ 1 )=SGN(Vd(L+l )-Vde)*(-1 . )+Vs<K) 

3130 

Vd ( Leo ) -Vs ( K+ 1 ) 

3140 

K-K+l 

3150 

GOTO  Loop2 

3160 

Increase : ! 

3161 

Vdcal=Vd(N+1 ) 

3170 

IF  ABS ( Vdca 1-Vde ) < . 0 1  THEN  GOTO  Out 

3180 

IF  SGN(Vdcal-Vde) >0  THEN  GOTO  Decrease 

3190 

Vsmin-Vs(K) 

3200 

GOTO  New  auess 

3210 

Decrease : ! 

3220 

Flg-1 

3230 

Vsmax-Vs(K) 

3240 

New  guess:  ! 

3250 

Vs(K  +  1  )  =  (Vsmax+Vsrriln)/2 

3260 

K  =  K  +  1 

3270 

Vd(Lco)=Vs(K) 

3280 

GOTO  Loop2 

3290 

Out:  ! 

3300 

RETURN 

3310 

! 

3320 

! 

3330 

i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  ii  i  i  i  i  i  i  i  i  i  m  M  i  m  i  i  i  i  i  i  i  i  i  i  i 

i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 

3340 

i 

j 

3350 

!  MODULE  D 

j 

3360 

j 

j 

3370 

!  CALCULATE  IB 

j 

3330 

1 

j 

3390 

■  i  i  i  i  i  i  i  i  i  i  ii  i  ii  i  i  i  i  i  i  i  i  i  i  ii  i  i  ii  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 

i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 

3400 

! 

3410 

Moduled:  ! 

3420 

Iocal-0 

3430 

FOR  L-l  TO  N 

3440 

IF  (C-Vd(L  +  D+Vbe)/(Kb*T))<-50  THEN  GOTO 

Around5 

3450 

Id -FN Inca ( Length ,  W  id th  .  N, L  )  +  <  Ies*(  1  -  At' )  * ( EXP (Q*Vbe/ ( Kb*T ) )  -1  )  +  ( Ics 

-At>IesH(EXP(0*(-Vd(L+1  )+Vbe>/<Kb*T))-1  ) ) 
3460  GOTO  3490 

3^70  AroundS:  ! 


J-3 


3480 

Ib=FN  Inca  (Length.  Width, N. !_)■*(  Ies*  (  i  -Af  )  *  ( EXP  ( Q*Vbe/  ( Kb-*T  )  1-1  ) - ( Ics 

-At  * 

Iesn 

3490 

Ibi(L)=Ib 

3500 

Ibcal=Ibcal+Ib 

3510 

NEXT  L 

3520 

RETURN 

3530 

i 

3540 

F  ini : 

i 

3550 

i 

3560 

'subend 

3570 

1 

3580 

1 

3590 

1 

3600 

DEF 

FNIncadength. Width. INTEGER  N. INTEGER  L) 

3610 

F=N 

3620 

Inca 

=  I/F 

3630 

RETURN  Inca 

3640 

FNEND 

3650 

i 

3660 

DEF 

FNWiddength, Width, INTEGER  N, INTEGER  L) 

3670 

Wid= 

Width 

3680 

RETURN  Hid 

3690 

FNEND 

3700 

i 

3710 

i 

rataf 


